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ABSTRACT 

We estimate the absolute magnitude distribution of galaxies which lie within a few 
hundred kiloparsccs of Mg II absorption systems. The absorption systems themselves 
lie along 1880 lines of sight to Sloan Digital Sky Survey Data Release 3 QSOs, have rest 
equivalent widths greater than 0.88A, and redshifts that range from 0.37 ^ z ^ 0.82. 
Our measurement is based on all galaxies identified by the SDSS photometric pipeline 
which lie within a projected distance of about 900 h^ 1 kpc of each QSO demonstrating 
absorption. The redshifts of these projected neighbors are not available, so we use a 
background subtraction technique to estimate the absolute magnitude distribution 
of true neighbors. Our method exploits the fact that, although we do not know the 
redshifts of the neighbors, we do know the redshift of the absorbers. The absolute 
magnitude distribution we find is well described by a bell-shaped curve peaking at 
rest-frame Mb ~ —20.2, corresponding to L/L* « 1.4. A comparison of this observed 
distribution to expected shapes derived from the literature suggests that it is unlikely 
to be drawn from a population that is dominated by late-type galaxies. However, when 
the sample is divided in half on the basis of rest-frame equivalent width, we find that 
some stronger equivalent width systems appear to be associated with later galaxy 
types. We use these shapes, along with the observed covering fraction of « 8%, to 
estimate the extent of the Mgll distribution around a galaxy. For an galaxy, this 
scale is about 70/i _1 kpc. We provide an analytic description of the method, which 
helps build intuition, and show that it is generally applicable to any datasct in which 
redshifts are only available for a small sub-sample. Hence, we expect it to aid in the 
analysis of galaxy scaling relations from photometric redshift datasets. 
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1 INTRODUCTION 

QSO absorption line systems have been the subject of nu- 
merous studies since their discovery and identification in 
the l a te 1960s jBachallll96Sl; i Burbidge, Lvnds. fc Stockton! 
1 19681 ; iBachall fc Spitzerl ll96S 1. Historically, these systems 
have been identified in spectra taken from the ground; at 
high enough redshift, atomic transitions with lines in the 
UV are redshifted into the atmospheric optical window, and 
indeed many such spectral lines have been used to identify 
these systems. Detailed studies of the number and kine- 
matics of these lines have greatly aided our understand- 
ing of the physical envir onment of the gas in these systems 
JChurchill fc Vogtll200lfl and have provided constraints on 
the amount of neutral gas in the univer se at high redshift 
|Prochaska. Herbert-Fort, fc Wolfel [20051 ) . The A2796A2803 
doublet of singly ionized magnesium (Mg II) is a popular 
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target of spectroscopic searches due to its relative ease of 
identific ation in spectra and its association with neutral hy- 
drogen (|Rao fc Turnshekll2000h . This makes it an excellent 
probe of neutral gas, particularly at redshifts below which 
Lyman alpha absorption is still outside the window observ- 
able from the ground. 

The connection between Mg II absorption sys- 
tems and luminous gal axies has been well established 
(|Bergeron fc BoisseHl99 J ), and models which place the ab- 
sorbing gas in the haloes of such galaxies have had some 
success in explaining the absorption characteristics seen 
jMo fc Miralda-Escudelll996l ; ISteidel et ai]|2002l ; iLin fc Zoul 
l200lf ). However, a more detailed connection between the ab- 
sorption systems and galaxy morphology as well as that 
between the absorption systems an d location with i n the 
galaxies is still uncer tain (but see iTinker fc Chenl 120081 : 
IChelouche et al.l 120081 . for more recent work). Deep imag- 
ing of 38 fields of QSOs that demonstrate Mg II ab- 
sorption in their spectra, combined with high resolution 



2 M. A. Caler, R. K. Sheth & B. Jain 



spectra of galaxies found in these fields, reveal the host 
galaxies to be mostly spiral galaxies, many with morpho- 
logical asymmetries suggesting a history of mild 



tational interactions ( Churchill. Kacprzak. fc Steide 
iKacprzak et all 120071 , 



gravi- 
1 2005; 



2008f ). Some fully saturated absorp 



tion systems have been shown to c orrespond to Damped Ly- 
man Alpha absorption systems (|Rao. Turnshek. fc Nestor] 
l2006h and hence to galaxies with a wide range of mor - 
phologies (jBowen. Tripp, fc Jenkinsll200ll ; iRao et al.ll2003h . 
whereas little is known about the host s of th e very weakest 
systems. ISteidel. Dickinson, fc Perssonl (| 19941 ) made the first 
measurement of the luminosity function of Mg II host galax- 
ies; their estimated K-band lu minosity function was found 
to be consistent with that of iMobasher. Sharpies, fc Ellisl 
l| 19931 ). with best-fit Schechter function parameters <f>* = 
3.0±0.7xl0" 2 {h/Mpcf,M* K = -25.1±0.3, a = -1.0±0.3 
for a sample of 58 galaxies. They also determined that the 
average absorber appears to be consistent with an Sb type 
galaxy (0.7 LJ), but noted a large spread (factor of ~ 70) 
in luminosity for the sample. 

Large surveys such as the Sloan Digital Sky Sur- 
vey (SDSS) have greatly increased the number of re- 
liably detected Mg II absorption systems. Searches 
for Mg II absorption systems within the spectra of 
SDSS QSOs have yielded close to ten thousand sys- 
tems for further study (iProchter. Prochaska. fc Buries! 



20061; 



Menard et al 



Nestor. Turnsek. fc Raol 
2008^ " 



2005; IBouche et al.l 



20061 ; 



However, progress in more detailed 
analyses of these systems is hindered by the shallowness of 
the photometry of these QSO fields and a lack of follow-up 
spectroscopy, as the SDSS is limited in its ability to take 
spectra of objects located within 50 arcseconds on the sky 
of each other. It is not feasible given current resources to 
carry out detailed follow up observations of thousands of 
fields; hence other methods must be used to gain informa- 
tion concerning the properties of the host galaxies of these 

systems and their environments. 

Recently IBouche, Murphy, fc Perouxl (|2004 l and 
IBouche et atl (|2006T l have described the use of cross- 
correlation techniques for studying the environments of 
Mgll absorbers; they estimate the host halos of absorbers 
to be ~ 5 x 1O 11 M0, and find an anti-correlation between 
a system's measured equivale nt width and the mass o f 
the halo of its host galaxy. IZibetti et all |2005l . 120071 ') 
have considered image stacking as a way to investigate 
the photometric properties of Mgll system host galaxies. 
Their stacking analysis provides an estimate of the average 
luminosities and colors of these galaxies. They conclude 
that the weaker absorbers are hosted by red, passively 
evolving galaxies, whereas the stronger absorbers are hosted 
by more actively star-forming galaxies. Our study is, in 
some ways, complementary to theirs. Both studies use 
SDSS photometric data to investigate the environments of 
Mgll systems These galaxies can then be used to determine 
the luminous environment of Mg II systems. Both of these 
methods provide a way to constrain properties of Mgll 
system host galaxies without any follow-up observations, 
though our method has the advantage of being somewhat 
easier to implement. 

In this work, we describe the results of an investigation 
into the absolute magnitude distribution of galaxies found 
near Mgll absorbers. Although we use the SDSS photomet- 
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Figure 1. Rcdshift distributions for our final sample of Mg II 
absorption systems and of the background QSOs. 



ric reductions to identify galaxies around QSOs with ab- 
sorption systems, these objects are generally too faint to 
have been part of the SDSS spectroscopic survey. Hence, al- 
though we have colors, we do not have redshifts for these 
galaxies. To compensate for the lack of redshift information, 
we use a background subtraction technique to isolate the 
galaxies physically associated with the Mg II host galaxies. 
We provide a discussion our sample in § [5] and describe our 
measurement technique in § [3] We present our results in § 21 
and summarize our findings and their implications in § [5] 
Throughout this paper we assume a ACDM cosmology with 
n M = 0.3, Q,x = 0.7 and that H = 100ft km s^ 1 Mpc" 1 . 



2 THE SAMPLE 
2.1 The absorbers 

Our sample of Mg II absorpti o n line systems comes from 
IProchter. Prochaska. fc Burled (|2006h . Full details of the 
sample selection method can be found there; we only give 
a brief summary here. Objects spectroscopically identified 
as QSOs in the SDSS Data Release 3 (DR3) are searched 
for evidence of Mg II absorption. The search is confined to 
QSOs with z > 0.35. A continuum fit for each spectrum is 
made using a 6-spline to fit the underlying QSO spectrum 
and a principal component analysis to fit any QSO emission 
lines. Spectral features are identified using a Gaussian fil- 
ter method; 3.5<r features are considered significant. Mg II 
lines are identified from the resulting list of lines by look- 
ing for lines matching the doublet separation. Features with 
measured equivalent width W / A2796 > 0.8A are compiled into 
the Mg II absorber sample. From searching a total of 46420 
QSOs in the SDSS DR3, there are a total of 9542 absorption 

systems in the final catalog. 

The absorbers in the IProchter. Prochaska. fc Burled 

|2006l ) sample span the equivalent width range 0.8 A to 5.0 
A, with a few detections out to 10 A. As Mg II absorption 
systems have been de tected to equivalent widths of 0.02 A 
(Churchill et al.lll999f ). we here are investigating the prop- 
erties of the strongest absorption systems. 

The redshift range over which these systems are de- 
tected is likewise broad; they are found over the full sensi- 
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tivity range of the SDSS spectrograph to the Mg II doublet 
lines, namely z = 0.35 to z = 2.2. However, the photometric 
catalog of the SDSS — on which we rely to study the galaxy 
neighbours of the Mg II host galaxies — is sensitive to galax- 
ies out to a redshift z ~ 1. For this reason, we choose to 
concentrate on the lowest redshift absorbers. We divide the 
full sample roughly into thirds and chose the lowest redshift 
bin for this study; we consider here absorbers in the redshift 
range z=0. 368-0. 820. A total of 2282 absorption systems fall 
within this redshift range. 

To ensure that we accurately investigate the environ- 
ment of the absorbing systems, we eliminate from our sam- 
ple all QSOs that show evidence for multiple systems in their 
spectra. This is necessary because we do not have redshift 
information for the majority of the galaxies in the neigh- 
borhood of the QSO position, and for lines of sight with 
multiple absorption systems it would be impossible to tell 
what galaxies were in the neighborhood of which absorber. 
This eliminates 142 systems from our sample, leaving 2140 
systems. We further eliminate all QSOs whose redshifts do 
not allow for possible detection of Mg II systems over the 
full redshift range z=0. 368-0. 82; this removes the lowest red- 
shift QSOs from our sample. This is done to eliminate possi- 
ble incompleteness effects in our absorption system sample. 
The redshift distributions for our absorption systems and 
the QSOs whose spectra they were found in is illustrated 
in figure [T] Our final sample is comprised of a total of 1880 
absorption systems. 

2.2 Reference sample 

The background subtraction technique (to be detailed in §[3]) 
we will use in this paper requires the construction of a sam- 
ple of random lines of sight to compare with the absorption 
systems. We construct this reference sample as follows. For 
each QSO whose spectrum demonstrates Mgll absorption 
(hereafter referred to as "absorbing QSOs"), three confirmed 
QSOs from the SDSS DR3 which do not demonstrate evi- 
dence for Mgll absorption along their line line of sight are 
chosen, each having similar redshifts and r— magnitudes to 
the absorbing QSO. These QSOs shall be hereafter referred 
to as "reference QSOs" . These are chosen from a set of 21543 
QSOs which have zqso > 0.82 and do not have a z > 0.36 
Mgll system along the line of sight. 

The redshift distributions of the absorbing and refer- 
ence QSO populations are shown in the top panel of figure[5] 
the r— magnitude distributions of the two samples are shown 
in the bottom panel. The matching of redshift and r magni- 
tude between the absorbing QSOs and reference QSOs en- 
sures similar S/N in their spectra, as well as similar spectral 
coverage; that is to say, the reference and absorbing QSOs 
have the same redshift window over which to detect Mgll, 
but the reference QSOs did not encounter an absorber. Each 
reference QSO is assigned a mock absorber whose properties 
are equal to those of the Mgll system found along the line 
of sight to the absorbing QSO for which it was selected to 
match. As this assigned system is a ghost, its properties will 
be uncorrelated with galaxies found in the imaged field of 
the reference QSO. 

The set of QSOs from which a reference sample can be 
drawn is about ten times larger than the sample of QSOs 
which contain absorbers - the ratio of lines of sight with 
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Figure 2. Redshift (top) and apparent m r magnitude (bottom) 
distributions of the absorbing and reference QSO populations. 



absorbers to those without is 1880/(1880 + 21543) = 0.08. 
If we include lines of sight with more than one absorption 
system, then this ratio is 0.085; we will return to this later. 



3 BACKGROUND SUBTRACTION 
TECHNIQUE 

Our method for estimating the absolute magnitude distri- 
bution of galaxies neighboring Mg II absorption systems 
closely follows the b ackground subtraction technique which 
lHansen et al.l (|2005l ) used when estimating galaxy cluster lu- 
minosity functions. We outline our technique here; a more 
detailed description of it can be found in the Appendix. 

We begin by using the SDSS (DR3) to find those ob- 
jects which are classified as galaxies and which lie within 
3 arcminutes of the 1880 QSOs of our Mg II absorption 
sample. We consider only galaxies with angular separations 
from the QSO position greater than 2 arcseconds, to elimi- 
nate any blending or seeing effects. Each galaxy is assigned 
the redshift of the absorption system associated with the 
QSO on which the field is centered. Angular separation is 
converted to comoving distance using a ACDM cosmology 
with Q.m = 0.3, S1a = 0.7. Due to the broad redshift range 
of our sample (Az = 0.45), 3 arcminutes corresponds to 
different comoving distances from the absorber host galaxy. 
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(For the mean redshift of our sample, z = 0.594, 2 arc- 
seconds corresponds to 14.8 h~ x kpc and 3 arcminutes to 
1.33 h^ 1 Mpc.) Higher redshift absorption systems sample 
galaxies to larger comoving separations than do lower red- 
shift systems; thus, our sample is incomplete at these large 
distances. We therefore consider only the subset of objects 
which lie within the range accessible over the entire redshift 
range: this fully sampled annulus spans comoving distances 
19.3 /i _1 kpc ^ d scp ^ 878 ft _1 kpc from the central QSO. 

We also use the redshift of the absorber to assign ab- 
solute magnitudes to each of the galaxies in its field. We 
set 

where cLl{z) is the luminosity distance to the galaxy and A 
is th e correction for extinction due to dust in the Milky Way 
from lSchlegel. Finkbeiner. fc Davis! (|l998r i. (Note that we do 
not include a k-correction term in our absolute magnitude 
calculations.) For those galaxies truly in the neighborhood 
of the Mg II absorption system, this procedure yields the 
true absolute magnitude. It of course yields an incorrect 
magnitude for all the other objects. 

We then follow the same procedure for each of the 
5640 reference QSO positions: Galaxies projected within 3 
arcminutes of each reference QSO are found and assigned 
redshifts as described above. Their angular separation from 
the reference QSO's position are converted to comoving dis- 
tances based on the redshift of the ghost absorber, and 
galaxies located within the fully sampled annulus are kept. 
These are assigned absolute magnitudes on the basis of their 
(ghost absorber) redshift. In this case, essentially all dis- 
tances and luminosities are incorrect. 

We now have absolute magnitude distributions centered 
on the absorber and reference populations. For the absorber 
population, 

iVabsorbcr (M) = ^neighbors (M ) + Arandom (M) J (2) 

whereas for the reference population, 

-^reference (-^0 ^random 

(M). (3) 

Here Aabsorber (M) denotes the number of galaxies with ab- 
solute magnitude M found in the field of an absorbing 
QSO, Af ran dom(A^) the number of such galaxies randomly 
projected into the field, and Aneighbors (M) the true neigh- 
bours of the absorption system host. If we subtract these 
two distributions — the absorber and reference distributions- 
taking care to account for the fact that we have three times 
as many QSOs in the reference catalog as in the absorber 
one, all that will remain is the contribution to the counts 
from those galaxies which are the true neighbors of the ab- 
sorption system. These are precisely those objects for which 
distances and absolute magnitudes were appropriately calcu- 
lated. Hence, the absolute magnitude distribution of Mgll 
system neighbour galaxies can be determined from this cal- 
culated difference in measured distributions. 

The Appendix provides a detailed calculation showing 
that the quantity which the background subtraction tech- 
nique actually estimates is 

-^neighbours (AT) ps V $ f dz abs £L 0(M|z abs ); (4) 



here </>(M|z a b s ) is the luminosity function in fields of effec- 
tive volume centered on absorbers (c.f. equation lA7|) . The 
Appendix also describes the results of testing our procedure 
on a mock catalog of galaxies, subjected to the same "ob- 
serving" limits as the SDSS, to ensure that it does in fact re- 
cover the absolute magnitude distribution of Mgll absorber 
neighbours. 

Our background subtraction technique relies on the ex- 
istence of an over density of galaxies around absorption 
systems which is not present in a random sample. Be- 
cause galaxies cluster over scales of ~ 100 kpc — 10,000 
kpc, such an over density of galaxies is expected. Hence, 
galaxy clustering ensures the viability of our technique. 
However, we have not yet addressed the question of how 
many lines of sight are required to make a statistically sig- 
nificant measurement. If the number of galaxies correlated 
with absorbers is some fraction C of the background counts, 
N Eibsorbcr (M) = (1 + C) iV ra ndom(Af). Then the 'signal' in 
our background subtraction method is C A ran d om (M), where 
Af ra ndom is proportional to the number of lines of sight. If 
the number of reference QSOs is many times that of ab- 
sorbing QSOs, or if the reference population was calculated 
analytically (c.f. equation lA5[l . then the Poisson noise on our 
measurement would be \/(l + C)A r ran dom(M). So if we want 
C AT ralldom (M) ^ 3 + C)A random (M) then the required 
sample size is 

N iandom {M) > 3 2 (1 ^ 2 C) . (5) 

This shows that our technique requires A r ran( jom(Af ) 3> 1 
if C < 1. For example, if C w 10%, then Abandon. ~ 10 3 . 
If there are only a few galaxies per line of sight, then our 
technique requires many lines of sight if C <C 1. If there 
are many galaxies per line of sight (e.g. if Mgll absorbers 
were centered on massive galaxy clusters), then fewer lines 
of sight would be required. 

Above, we noted that we did not include a k-correction 
term when calculating absolute magnitudes for galaxies in 
our sample. As a result, our observed r— magnitudes do not 
correspond to the rest-frame r— magnitudes of Mgll system 
neighbour galaxies. For most of the absorption systems in 
our sample, the r— band absolute magnitudes we calculate 
are close to rest-frame B-band absolute magnitudes; systems 
with z > 0.7 are closer to rest-frame U. The SDSS is not 
expected to find many galaxies with z > 0.7, so we do not 
expect to find many Mgll system neighbour galaxies in this 
redshift range. Thus including them in our survey should 
not contaminate our results very much. 

Given the background subtracted absolute magnitude 
distribution, we could calculate a luminosity function once 
we account for the fact that our galaxy survey is magni- 
tude limited, so more luminous galaxies are seen to greater 
distances. This could be done by carefully implementing 
Schmidt's V ma x method for the particulars of our survey (see, 
e.g., Appendix). In practice, because our sample spans a rea- 
sonable range in z, and we do not apply k-corrections, we 
do not really distinguish k-corrections from evolution; there- 
fore, estimating the luminosity function using such a Knax 
method is slightly tricky. For this reason, we present our 
results in terms of the observed absolute magnitude distri- 
bution, -/V nc i g hbours(Af), rather than an evolution corrected 
luminosity function. 
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Figure 3. Absolute magnitude distributions for the absorber 
population, plotted in black, and the reference population, plot- 
ted in blue. 
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Figure 4. Difference between the reference samples. A line at 
zero counts is shown for reference. 



4 RESULTS 

Before carrying out the background subtraction procedure 
detailed in §[3] we first present the absolute magnitude dis- 
tributions of the absorber and reference populations in fig- 
ure [3] to determine if there is a non- negligible difference be- 
tween the two. Such a difference would correspond to the 
true neighbours of the Mgll system host; this is the pop- 
ulation we aim to isolate. Figure [3] shows that there is an 
excess of absorber population counts over reference popu- 
lation counts from —22.0 ^ M r ^ —18.0. This indicates, 
even before our method has been applied, that we have suc- 
cessfully found an over density of galaxies in the absorber 
population, not extant in the reference population, which are 
associated with the absorption systems in their fields. These 
galaxies are all that should remain upon implementation of 
our background subtraction technique. 



4.1 Reference sample 

We test that our reference population is truly comprised of 
galaxies randomly projected into the fields of the central ref- 
erence QSOs before implementing our technique. To do so, 
we split our reference catalog into three equal-sized catalogs, 
each containing one reference QSO for every absorbing QSO. 
We then subtract the mutual absolute magnitude counts for 
each catalog pair to see if the resulting distribution is cen- 
tered on zero to within the scatter. Figure [3] shows the vari- 
ous pairwise differences between the three reference catalogs. 
The absence of any real feature in this figure is reassuring; 
it suggests that our reference sample is truly random. 

We perform one additional test on our reference sam- 
ple to ensure its random nature. As described in the Ap- 
pendix, the expected distribution of apparent magnitudes 
in a randomly chosen field can be computed analytically if 
the luminosity function is known (equation IA1|) . The result 
of assuming that this is given by the COMBO-17 survey 
i|Wolf et allbOOSt ) and the ESP survey (IZucca et al.)ll997l ) 



(converted to our cosmology following iLiske et al.l ((20031 )1 
is shown by the dashed curve in the top panel of Figure [S] 
(These surveys span the same redshift range as our dataset.) 
The histogram shows the actual counts; while our calculated 



apparent magnitude distribution matches the observed one 
through m r ~ 21.5, we over-estimate the number of galax- 
ies with m r > 21.5. This is due, in part, to our neglect of 
k-corrections, which have the effect of smearing out what 
would otherwise be a sharp cut at m r . The reference sam- 
ple for each absorber is constructed by shifting this appar- 
ent magnitude distribution by factor which depends on z a b s - 
Since the distribution of z a b s is known, we can convert this 
apparent magnitude distribution to Af re f erence (M) (see equa- 
tion \K4§. The dashed line in the bottom panel of Figure [5] 
shows the result of transforming the dashed curve in the up- 
per panel in this way (i.e., using equation I A4|l . and the his- 
togram shows the actual distribution of Af re ference(M). The 
small differences between the predicted and actual distri- 
butions can be traced to the differences in the top panel. 
Given our understanding of the nature of these differences, 
and the overall agreement in shape between the calculated 
and observed distributions, we conclude that our reference 
population is indeed consistent with a population of galaxies 
randomly projected into a field. 



4.2 Background subtraction and the full sample 

We now implement the procedure outlined in § [3] to 
isolate Mgll system neighbour galaxies. Figure [S] shows 

cf'crcncc 

(M) as a function of M, for each 
of the three equal-sized reference catalogs we used to test 
our reference sample (see Figure |3| . In all cases, there is 
a broad well-defined bump centered on M r ~ —20.2. That 
we see a bump of approximately the same shape peaking 
at approximately the same place in all three distributions 
indicates that we measure a real underlying one and not a 
by-product of noise. Moreover, the bell-like shape is simi- 
lar to that which one typically obtains if the raw galaxy 
counts in an apparent magnitude limited survey are plotted 
as a function of absolute magnitude. This gives us confidence 
that our method has worked. In what follows, we will reduce 
the noise in our estimate of the background counts by using 
the full reference catalog compiled in §[3] when subtracting, 
rather than the three equal-sized catalogs. When we do so, 
we detect a total of 2797 galaxies — about 1.5 times the 
number of absorbers. 
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Figure 5. Apparent (top) and absolute (bottom) magnitude dis- 
tributions (histograms). The dashed lines show predictions based 
on equations (I AH and \A4\ for the top and bottom panels re- 
spectively. 



Because of luminosity evolution in our galaxy popula- 
tion over the range 0.37 ^ z ^ 0.82, further complicated 
by potential confusion with k-corrections (which we ignore) , 
we do not convert our background subtracted counts into 
an estimate the luminosity function and show the resulting 
(f>(M). Rather, we compare these counts directly with pre- 
dictions based on inserting various models for the luminosity 
function and its evolution into equation Q. When doing so, 
we treat as a free parameter; recall that this parameter is 
a measure of the effective volume associated with absorbers. 

As we did when modelling the reference sample, 
we use luminosity functions from the COMBO-17 sur- 
vey to model our background subtracted sample. Since we 
do not k-correct our absolute magnitudes, our observed 
r— magnitudes correspond to their rest-frame B-magnitudes 
up to z ~ 0.7, covering almost the entire redshift range of 
galaxies considered (as noted above). The redshift range of 
this survey is well-matched to that of the galaxies we detect 
with our method; galaxy luminosity functions for four dif- 
ferent populations are reported therein for mean redshifts 
z=0.3, z=0.5, z=0.7, z=0.9, and z=l.l. This amply covers 
the redshift range over which we detect Mgll neighbor galax- 
ies. 

As noted above, COMBO-17 reports luminosity func- 



Figure 6. Difference in observed counts between the absorber 
sample and each of the three reference samples. 
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Figure 7. Background subtracted absolute magnitude distribu- 
tion of Mgll system neighbour galaxies (histogram) and expecta- 
tions based on inserting luminosity functions from the COMBO- 
17 survey into equation flt . 



tions for four galaxy types, as well as an estimate for the full 
population. In our estimate of the reference counts, we used 
only the full population. Here, to explore the possibility that 
Mgll absorption system neighbors may be preferentially of 
one type, we show the result of inserting each of the four 
luminosity functions, as well as that for the full population, 
into equation (j4|. The four smooth curves in Figure [7] show 
these different predictions; in each case the overall ampli- 
tude has been left as a free parameter. The histogram 
shows the background subtracted counts that we measure. 
Evidently, the COMBO-17 Type 1 luminosity function pro- 
vides the best description of our measurements. Both the 
Type 2 and Type 3 luminosity functions provide an accept- 
able fit to the faint-M end of the distribution, but neither 
match well at the bright end. The same is true for the total 
(i.e. all types) luminosity function. The Type 4 luminosity 
function provides the poorest fit to the data. 

We conclude that our observed distribution seems to be 
most consistent with that calculated for a COMBO-1 7 Type 
1 luminosity function, which according to I Wolf et all (|2003T ) 
is drawn from galaxy types ranging from E — Sa. Taken at 
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Figure 8. Absolute magnitude distributions for the weak 
and strong (bottom) absorber and reference populations. 



(top) 



face value, then, our comparison suggests that the neigh- 
bours of Mgll absorption systems are unlikely to be very 
late-type galaxies. It is interesting to note, in passing, that 
we find that we require as 355 to fit our background 
subtracted counts when the Type 1 luminosity function is 
considered. This in equation (|A11|) suggests a correlation 
length ro of about 4.6 /i -1 Mpc, which is not unreasonable. 



4.3 Samples split by equivalent width 

To investigate the possibility that absorption systems of dif- 
ferent strengths are associated with different environments, 
we divide our sample in half according to equivalent width. 
The dividing point occurs at rest-frame equivalent width 
(REW) 1.28A. Hereafter, we shall refer to the sub-sample 
with 0.88 REW sC 1.28 A as the "weak" sample, and 
the sub-sample with REW > 1.28 A as the "strong" sam- 
ple. Note that, when we divide the sample, each absorber 
QSO "keeps" its three reference QSOs; this ensures that the 
redshift and apparent magnitude distributions of our refer- 
ence samples match those of the subsamples. 

As we did with the full population (compare Figure [3]) , 
we first present the absolute magnitude distributions of the 
weak and strong absorber and reference populations. Both 
the top and bottom panels of Figure [S] show an excess of 
absorber population counts over reference population counts 
over the range —22.0 sC M r —19.0. This indicates that we 



Figure 9. Background subtracted absolute magnitude distribu- 
tions of the weak and strong sub-samples. The background sub- 
tracted absolute magnitude distribution of the full sample is also 
shown for reference. 



find neighbours around both weak and strong absorption 
systems rather than around just one sub-sample. 

The light-grey and dark-grey histograms in Figure [9] 
show the result of subtracting reference from absorber 
counts in the weak and strong samples, respectively. Note 
that, because of how we have defined our reference samples 
(i.e., each absorber QSO "keeps" its three reference QSOs), 
the background subtracted counts sum to give the absolute 
magnitude distribution of the total sample (black histogram, 
same as shown in Figure^, by definition. Both distributions 
are clearly different from zero; they both show a broad bump 
over the absolute magnitude range —19.0 M r sj —22.0. 
Both the weak and strong distributions appear to peak at 
M r m —20.2, though the peak of the strong distribution is 
slightly better defined than the peak of the weak distribu- 
tion. 

The two distributions do, however, appear to be slightly 
offset; the strong distribution extends from —18.0 M r sj 
—22.0, whereas the weak distribution extends from —19.0 si 
M r sj —23.5. The weak distribution also features a high- 
M tail which is not present in the strong one. As another 
indication of the difference in distributions, note that the 
background subtracted absolute magnitude distribution for 
the weak sample contains 1679 galaxies, whereas that for 
the strong sample contains 1118 galaxies - even though 
the two samples were centered on the same number of ab- 
sorbers: 940. This suggests that galaxies near weak systems 
are more readily detectable by the SDSS than are galax- 
ies near stronger systems, which would imply that weaker 
systems are surrounded by brighter galaxies provided both 
sub-samples have similar redshift distributions. 

If the redshift distributions of the two absorber popu- 
lations were the same, we could conclude from the absolute 
magnitude distributions shown in Figure[9]that the luminos- 
ity function of galaxies centered on strong absorbers includes 
more faint objects. Figure [10] suggests that dN/dz is indeed 
similar for the two populations. As an extra check, we per- 
formed the same analysis as before - using equation (|A7[l to 
calculate absolute magnitude distributions from COMBO- 
17 luminosity functions, using the redshift distribution ap- 
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Figure 10. Rodshift distribution of weak (red) and strong (blue) 
Mgll absorption systems. 



propriate for each sub-sample, and leaving Vj as a free pa- 
rameter. 

The top panel in Figure [TT] shows results for the weak 
absorber population; a COMBO-17 Type 1 luminosity func- 
tion provides a good description of our measurements, at 
least over the range —22 ^ M ^ —18. The bottom panel 
shows that although a Type 1 luminosity function also pro- 
vides the best description of the bright end, there is a no- 
ticeable excess at low luminosities, which could be described 
by a Type 4 luminosity function. Thus, it appears that 
whereas weak absorbers tend to not be associated with late- 
type galaxies, at least some strong ab sorbers may well be . 
Thi s is consistent with the results of IZibetti et al.l (|2007l ) 
and lBouche et al.l (|2006t ). who used very different techniques 
from ours. 




stronq b.s. counts 
:':-17 T' LF 
' I BO-17 T2 LF 
C0MB0-1 7 T3 LF 
BO-17 T4 LF 
•V'\ BO-17 ALL LF 




Figure 11. Background subtracted absolute magnitude distribu- 
tions for the weak (top) and strong (bottom) sub-samples, com- 
pared with the expected distribution based on various luminosity 
functions from the COMBO-17 survey. 



5 DISCUSSION 

5.1 Summary of results 

We have estimated the absolute magnitude distribution 
of galaxies which lie within about 1 /i _1 Mpc of Mgll 
absorption line systems. Our sample of 1880 absorbers, 
which is drawn from the SPSS D R3 Mgll catalog of 
IProchter. Prochaska. fc Burled (|2006h . spans the redshift 
range 0.368 ^ z ^ 0.820, and consists of systems with 
rest-frame equivalent width REW > 0.88A. Lines of sight 
demonstrating Mgll absorption at multiple redshifts have 
been eliminated, as have QSOs whose spectra would not al- 
low systems to be detected over the entire redshift range. 

Most of the galaxies in SDSS imaging which lie within 
about 1 7i _ Mpc of each of these objects have five band 
imaging but no spectra. Hence their redshifts are not known; 
this, of course, complicates estimates of the absolute mag- 
nitude distribution. In principle, w e could estim ate lumi- 
nosities using photometric redshifts; IShetbl l|2007l ) describes 
how to derive accurate estimates of the luminosity function 
from photo-zs. In the present context, however, we use only 
the photometric information for the galaxies, because we 
do know the redshift of the absorber. Hence, we use this 
information and a background subtraction technique to sta- 
tistically remove foreground and background galaxies. 



Galaxies located within about 3 arcminutes of absorb- 
ing QSO positions in the SDSS DR3 imaging are assigned 
the redshift of the absorption system; their absolute magni- 
tudes are calculated based on this redshift. To correct for 
foreground and background galaxies that have been pro- 
jected into the field, we carry out a similar procedure on 
a reference set of QSOs. The reference QSOs have the same 
redshift and r-band magnitudes as the QSOs with Mgll ab- 
sorption features (Figure [5J), but do not demonstrate ab- 
sorption in their spectra. Galaxies in the fields surrounding 
a given reference QSOs are assigned the same redshift as the 
QSO with absorption features which that reference QSO was 
chosen to represent. A detailed analytic description of the 
method is provided in Appendix Al. Tests of our technique 
on a mock catalog of galaxies that is subjected to a similar 
reduction process as was used for the actual data indicate 
that it works well (Appendix |A)| . 

The background-subtracted absolute magnitude distri- 
bution we observe peaks at M r ~ —20.2 and is roughly bell- 
shaped over the range —22.0 ^ M r ^- —18.0 (Figure [6]). 
Comparison with expectations from published luminosity 
functions over this same redshift range (from COMBO-17) 
suggests that the galaxies within ~ 1 /i -1 Mpc of fields cen- 
tered on absorbers are of Type 1, meaning they may be 
ellipticals or Sa's (Figure [7}. Subdivision into weak and 
strong systems suggests that late-type galaxies may be asso- 
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ciated with the stronger systems (Figure [TT} , but are most 
likely not associated with the weaker systems, in agree- 
ment with previous work based on ver y different methods 
jBouche et alj|200d : lzibetti et al.ll2007l ). 



5.2 Implications for covering fraction 

The fraction of eligible lines of sight which contain an ab- 
sorber is about 8 percent. The galaxies which our back- 
ground subtraction technique identifies are the neighbours of 
absorption systems. If we assume that at least some of these 
neighbours are in fact the Mgll system hosts themselves, 
then this observed fraction can be used to place interesting 
constraints on the sizes of absorbers as follows. (Note that 
this is not the main focus of our paper). If absorbers are 
associated with galaxies, and are sufficiently rare such that 
they do not overlap, then the area on the sky they cover is 



(6) 



where <f>(L\z) is the luminosity function, cLa{z) is the an- 
gular diameter distance, R(L,z) is the radius out to which 
Mgll absorption is seen in the observed range of equivalent 
widths (e.g., some models have equivalent width decreasing 
with distance from the center of the galaxy, with a nor- 
malization which depends on L), and re is the fraction of 
galaxies which have an absorber (e.g., if the absorbers are 
small clouds within a galaxy, or if only a fraction of galaxies 
actually contain Mgll). 

It is common to parametrize 



R(L,z) = R„i 



(7) 



and to assume re is a constant ISteidel. Dickinson, fc Persson 
(1 19941) . Recent work s uggests that f3 « 0.35 ( Chen fc Tinker 
120081 ). or (3 « 0.2 (|Kacprzak et al.l 120081 ). Then, for a 
Schechter luminosity function with faint-end slope a, the 
expected covering fraction for our sample is: 



F = K(c/H o )(TrRi w )r(l + a + 20) 



I 

Jo 



dz 



0.37 N /n (l + z) 3 + (l-O„) 



(8) 



If the galaxies isolated by our background subtraction 
method are drawn from a COMBO-17 Type 1 z = 0.5 lu- 
minosity function (cf> t = 0.0028 ft 3 Mpc -3 and a = 0.52) and 
we ignore evolution, then the expected covering fraction is 

R*w 



F = 0.25 k 



100ft -1 kpc 



(9) 



We have used (3 = 0.35 in our calculations. Setting f3 = 0.2 
decreases the right hand side by ten percent. Since we 
find F = 0.08, re = 0.32 if R,w = 100ft -1 kpc, and 
R*w ~ 57ft -1 kpc if re — 1. We can account for luminos- 
ity evolution in our calculation of equation (|8]) by using the 
Type 1 luminosity functions reported by Wolf et al. (2003) 
for redshifts z=0.3, z=0.5, z=0.7, and assuming no evolution 
between 0.3 «S z < 0.5, 0.5 < z < 0.7, and 0.7 < z < 0.9. 
Doing so, we find that 



F = 0.18k 



R* 



100ft- 1 kpc 



(10) 



For F = 0.08 and R*w = 100ft -1 kpc, re = 0.44 is required; 
if instead re = 1 then R*w ~ 67ft" kpc. The values of re we 
obtain with our rough estimates are lower than those found 
in the literature, but not unreasonably so. Similarly, our 
estimate of R,w is reasonable. We note, however, that for 
a galaxy with REW ^ 0.88, a size of « 70ft -1 kpc is larger 
than would be expected from Figure 4 of IChen fc Tinker! 
(|200Sl ). 

To illustrate that the estimated k and R*w depend on 
the adopted luminosity function, suppose that these galaxies 
are drawn from the COMBO-17 total luminosity function 



0.018 ft 3 Mpc" 3 and a = -1.1 at z = 0.5). Then 



F = 2.05 k 



R*\ 



100ft- 1 kpc 



(11) 



if we ignore evolution. Setting (3 = 0.2 approximately dou- 
bles the right hand side. With F = 0.08, R*w = lOO/i -1 kpc 
implies re = 0.039, and R,w ~ 19.75ft -1 kpc if re = 1. Ac- 
counting for evolution as before (in this case a also evolves, 
so we keep the Gamma function piece inside the redshift 
integral) we find 

- I > (12) 



100ft -1 kpc 

implying re = 0.026 if R,w ~ 100ft -1 kpc and R,w ~ 
16.25ft -1 kpc if re = 1. These values are illustrative only, 
because this luminosity function does not result in good 
agreement with our background subtracted counts. 



5.3 Future work 

Our background technique is generally applicable to other 
studies in which redshifts are known for only a small sub- 
set of objects, but these are known to be correlated with a 
larger sample, for which we wish to estimate luminosities. 
For this reason, we provide a detailed analytic description 
of the method in Appendix Al. The analytic arguments al- 
lowed us to check a number of intermediate steps in our 
method (Figures [5l and lAlj) . and to estimate the required 
sample size for implementing this technique (equation 

It is also straightforward to extend our technique to 
study the clustering of the absorbers (whose redshifts are 
known) with galaxies in the SDSS photometric sample 
(whose redshifts are not known) to study the Mgll absorber- 
galaxy cross correlation function. This is the subject of work 
in progress. 
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APPENDIX A: MOCK CATALOG 
CONSTRUCTION AND TESTS 

To ensure that the background subtraction method we use 
in our luminosity function estimation yields robust, accu- 
rate results, we test it on a mock catalog of galaxies. This 
mock catalog will need to incorporate a model of galaxy 
clustering, for it is the clustering of galaxies over scales of 
« 100 kpc — 10,000 kpc which ensures the viability of our 
technique. Because galaxies cluster, there will be an over 
density of them, compared to random, which surround a 
selected galaxy and share its redshift, in addition to unasso- 
ciated field galaxies randomly projected into the field. This 
over density of galaxies is all that remains when we subtract 
away the random field galaxies. Were we to ignore galaxy 
clustering in our mock catalogs, our method would yield 
nothing, for all galaxies in the field of a selected one would 
be random, unassociated ones, and would thus be eliminated 
from our sample upon background subtraction. 

Due to the nature of our technique, we actually need 
to produce two mock catalogs: one catalog of mock galax- 
ies projected near simulated Mg II absorption line systems, 
and one of mock galaxies projected near simulated reference 
QSOs. While the same mock galaxy generating code will be 
used to construct both catalogs, the method by which simu- 
lated Mg II absorption line systems and simulated reference 
points are chosen must necessarily be different. In what fol- 
lows, we first present an analytical overview of our technique 
to give the reader greater intuition for how it works. We next 
describe our mock galaxy generation code, then detail how 
the simulated Mg II absorption systems and references, re- 
spectively, are chosen. Finally, we explain how both catalogs 
of simulated galaxies are compiled and present the results 
of applying our background subtraction technique to them. 
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Al Analytic calculation of background 
subtraction technique 

We begin our discussion with an analytic calculation which 
shows how our technique works. 

In a flux limited sample which covers some fraction / s k y 
of the sky, the observed number of objects with apparent 
magnitude m is 

N{m) = /sky^ dz -^j~T~ j dMcf>{M\z) 



X St 



D m 



M + 5 log + k(z) 



I sky 



dz 



dV(z) 
dz 



m — 5 lof 



d L («) 



10 pc 



- fc(z)Lj, 



where cj){M.\z 



10 pc 

is the luminosity function at z, dh(z) 



(Al) 



the luminosity distance to an object at z, and k(z) is its 
fc— correction. The surface density (number per unit area) 
of objects is 



n(m) 



N(m) 



(A2) 



47T /sky 

If we assign all these objects the same redshift (and k- 
correction), then (|A1[) will also describe the shape of the 
"luminosity" distribution which results, except for a con- 
stant shift. If we do this assignment for a number of differ- 
ent choices of redshift, the distribution of "luminosities" will 
be given by simply shifting this shape for each redshift and 
summing up the result. Thus, 

dN f m ™ 

■w(2abs) / dmn(m) 



N ran (M) = / dz ab . 



x Sr> (m = m - 



Gfc(Zabs) 



fc(2abe))- (A3) 



10 pc 

Here dN/dz^bs is the distribution of redshifts to be assigned 
to the objects (in this case the distribution of absorber red- 
shifts), and we have allowed for the possibility that angular 
sizes uj of fields associated with redshift z a bs ma y depend on 
z a bs- If we explicitly include in (IA3|) our expression (|A1[) for 
the surface density of objects, we obtain 

ax (M) 



Wran(Af) 



dN W(-Zabs) 



x 4> - 5 log 



djZabs 47T 

dh(z) 



dz 



JM) 



dV 

dz 



k(z) + fc(Zabs) 



(A4) 



dL(-Zabs) 

In essence, (|A4fl describes the expected distribution of ab- 
solute magnitudes of objects which have random angular 
positions in a field but the same redshift distribution as the 
absorber catalog. This is precisely what we want our catalog 
of mock galaxies projected near simulated reference QSOs 
to contain. We can simplify this expression further if, as in 
the case of our mock catalog, the luminosity function does 
not evolve and there are no fc-corrections: 



JVran(M) = / dz 



dN 

dz, hs 



(M) 

d L (z) 



(M > dV{z) 
dz — 
dz 



x^(A/-51og-^U 



(A5) 



We shall show in 



dL(Zabs 

that ()A5|) provides an excellent 



description of the absolute magnitude counts of simulated 
galaxies in our mock catalog counterpart to the reference 
population of § 

The above considers the case in which fields are centered 
on a random point on the sky. If instead fields are centered 
on objects which are correlated with other objects in the 
field, there will be an additional contribution which comes 
from this spatial correlation £. Namely, 



N e {M) 



Zmax(A/) dN 

dZabs j <^(M|z a bs 

Zmi„(M) 



X 



P r max f c 

2tv / dr p r p / 

J J — ( 



*/£(»>, J/|z a bs),(A6) 



where we assume that luminosity distances di,(z) and k- 
corrections do not change appreciably over the range of 
scales on which £ is not negligible, and that £ does not de- 
pend on luminosity. It is this extra term which the back- 
ground subtraction technique isolates. 

To gain intuition about this term, suppose that £ does 
not evolve over the redshift range spanned by the absorbers. 
Then the term on the second line of ()A6[) above is simply 
a constant (and, under the current hypothesis, independent 
of M). For example, if the correlation function had the form 
£(r) = {r /r) 2 , then 



N t (M) « V t 
where 



ax(M) ^ 

), (A7) 



V ( = 2n (7rrg) rmax rmin . 

ro 



(A8) 



Further, if we assume the luminosity function does not 
evolve over the range of redshifts spanned by z a bs, we could 
further reduce (|A7|) to 

/•Zmax(M) 

JV £ ( Af ) « V £ 0(M) / Cteabs (dJV/dZab B ) ■ ( A9) 



The right hand side of ()A9|) is proportional to the luminosity 
function times the number of fields lying in the redshift range 
wherein an object of absolute magnitude M would have been 
observed in a flux limited catalog. 

Let A^abs denote the total number of absorbers in a cat- 
alog. If we define 

/•Zmax(M) 

F abs (M) = N~£ / dz abs (dN/dz abs ), (A10) 

iz min (M) 



then 



N ( (M) » Vt 4>{M) N aha Fa bs (M). 



(All) 



For n abs = 0.001/(/i _3 Mpc 3 ), r = 5ft _1 Mpc, r max = 
l/i _1 Mpc and r min = 0.01/i _1 Mpc, we find that A 5 (M) w 
0.5 <f>(M) N^s F & b s (M). 

Notice that if the absorbers were uniformly distributed 
in comoving volume, i.e., dN/dz ahs = n a bs /sk y dV/dz^a, 
then 

N ( (M) n fiabs V e 4>{M) / sky [V max (M) - V min (M)]. (A12) 

If one then weights objects with luminosity M by the inverse 
of /sky [Knax(M) — Vmin(M)], then the resulting distribution 
will be proportional to the luminosity function 4>(M). The 
constant of proportionality is the product of the number 
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density of absorbers and the correlated volume. In general, 
the absorbers will not be uniform in comoving volume, e.g., 
because of the QSO redshift distribution, or S/N issues with 
the spectrograph. If so, the appropriate weighting factor is 
F ahB (M) rather than V m ^{M). 

A2 Mock catalog construction 

Having provided an analytical description of our technique, 
we are ready to detail the construction of our mock catalogs, 
beginning with our mock galaxy generating code. Before we 
doing so, however, we must consider how in our construction 
to place galaxies near each other so as to mimic a clustering 
signal. Ideally, this would be done by populating dark matter 
halos produced in N-body simulations with galaxies accord- 
ing to the halo distribution function. Since we do not have 
such simulations to work with, we chose to mimic galaxy 
clustering by placing galaxies into groups of fixed comoving 
spherical volume, each group containing 20 galaxies. This 
is not meant to be an accurate description of real galaxy 
clustering properties, but rather a toy model through which 
we can illustrate the effects of clustering. Within a group, 
galaxies are uniformly distributed within a spherical volume 
comoving radius 288/i _1 kpc (i.e., the comoving volume of 
each group is 0.1ft~ 3 Mpc 3 ). Projected onto the sky, the 
angular size of such a group at a redshift z — 0.82 is 30 
arcseconds; at a redshift z = 0.37 it is 1 arcminute, and at a 
redshift z — 0.1 it is 3.4 arcminutes. (The background cos- 
mology is ACDM with parameters given in the main text.) 

Having chosen a method for placing galaxies into 
groups, we are ready to populate a simulated volume with 
them. A 1.56 square degree patch of sky spanning the red- 
shift range ^ z ^ 1 forms the volume which we populate 
with mock galaxies. Galaxy groups are placed at random 
within this volume with number density 8.15 x 10~~ i h 3 Mpc 3 
and their redshifts calculated accordingly. Each group is 
then assigned a random position within this patch of the sky. 
The angular size on the sky of each group is computed (this 
depends on its assigned redshift); its 20 member galaxies are 
then each assigned its redshift as well as a sky position. Lu- 
minosities for e ach ga laxy are then drawn according to the 
iBlanton et al.l (|2003bT l luminosity function, evolved to the 
mean redshift (z = 0.594) of our M g II absorption sample by 
the prescription of lLin et al.l(| 19991 ). This evolved luminosity 
function has <j>* = 1.54 x 10 -2 (/i^Mpc)" 3 , M* = -21.44, 
and a = —1.04. We shall hereafter refer to this luminos- 
ity function as the "input luminosity function." Once all 
galaxies have a redshift and absolute r— magnitude, their 
apparent r— magnitudes can be determined. 

Each mock galaxy now has a redshift, an angular posi- 
tion on the sky, an absolute magnitude M r and an apparent 
magnitude m r ; this is all the data we need to select either a 
simulated Mg II absorption system or a simulated reference 
QSO and find all galaxies projected near it. We refer to this 
set of objects as the "full" mock catalog. A subset of these 
objects will satisfy the SDSS magnitude limit of m r — 22.5; 
we call this the "apparent magnitude limited" mock catalog. 

We must now choose a galaxy from the "full" mock cat- 
alog to serve as our mock Mg II absorber (selecting absorbers 
from the "apparent magnitude limited" catalog is inappro- 
priate, since detection of a galaxy by Mg II absorption is 
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Figure Al. Absolute magnitude counts per bin obtained from 
the procedure outlined in § [XT] and the mock reference neighbour 
catalog. The black histogram shows the counts from the mock 
reference neighbour catalog; the smooth red curve is calculated 
as in § IA1I The top panel displays the counts on a linear scale; in 
the bottom one, they are shown on a log scale. 



independent of its apparent brightness). We would like the 
set of such galaxies to match as closely as possible the red- 
shift distribution of the Mg II absorption systems in our 
data. To achieve this, we randomly select an absorber from 
the SDSS data set and choose a mock absorber from our 
"full" catalog that has the same redshift. We add one addi- 
tional requirement, which is that the distance to the nearest 
object in the "apparent magnitude limited" mock catalog 
match that of the nearest neighbour seen in the SDSS data, 
and that it has a similar apparent magnitude. Once our mock 
Mg II absorber has been chosen, we find those galaxies in 
the "apparent magnitude limited" catalog located within 3 
arcminutes of its position. Once found, they are re-assigned 
the redshift of the mock absorber and compiled into a final 
catalog, hereafter referred to as the "mock Mg II neighbour 
catalog". We repeat the process of finding a mock Mg II 
absorber and its detectable neighbours-that is, those neigh- 
bouring galaxies with m r ^ 22.5-until there is one mock Mg 
II system for each real one in our data sample, and all their 
neighbouring "detectable" galaxies have been added to the 
mock Mg II neighbour catalog. 

The next step is to create a counterpart to the reference 
sample. For each absorber with known redshift, we now se- 
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Figure A2. Absolute magnitude counts per bin obtained from 
the background subtraction of our mock Mg II neighbour and 
mock reference neighbour catalogs. The construction of these two 
mock catalogs is outlined in the text. The solid line is the expected 
shape of the distribution and is calculated as in § [X] 

lect a random sky position to serve as a mock "reference 
QSO". (A random sky position is chosen because our mock 
"reference QSO" position need not correspond to a mock 
galaxy's, as it did for the mock Mg II systems.) We en- 
sure that this random point does not coincide with that of 
a mock galaxy in the redshift range 0.37 < z < 0.82 of the 
"full" mock catalog, for in this case it could potentially have 
been a mock absorber. We then search for mock galaxies in 
the "apparent magnitude limited" subset which are located 
within 3 arcminutes of this position. These objects are all 
assigned the same redshift (that of the absorber), and then 
compiled into what we call the "mock reference neighbour 
catalog" . We repeat this three times for each absorber in our 
data set. 

To test that all has worked well, the histogram in Fig- 
ure [AT] shows the counts per absolute magnitude bin in the 
mock reference neighbour catalog described above. The solid 
red line shows our analytic calculation from the previous 
section (see IA5[) . The correspondence between the two is 
excellent. 

With our mock Mg II neighbour and mock reference 
neighbour catalogs now compiled, we apply the procedure 
outlined in § [3] to estimate the absolute magnitude distri- 
bution of the galaxies they contain. The result of doing so 
is shown in Figure IA2I The solid line in this figure is our 
analytic calculation from equation (|A9|) ). Though noisy, the 
shape of the absolute magnitude distribution mirrors the ex- 
pected one. This demonstrates that our method can success- 
fully recover the underlying absolute magnitude distribution 
of galaxies correlated with Mgll absorption systems. 

In the mock catalog, there are no fc-corrections and 
the luminosity function does not evolve. In this case, equa- 
tion (|A11[) suggests that the N(M) distribution can be 
turned into an estimate of cj>(M) simply by applying an ap- 
propriate M dependent weight to each galaxy. In practice, 
neither fc-corrections nor evolution can be ignored, so this 
complicates the relation between N(M) and <j>(M\z). This 
is why, in the main text, we present plots of N(M) but not 

of <t>. 



